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The  Conversion  of  Off-Vertical  Observations 
of  Total  Electron  Content  Into  Equivalent 
Vertical-Incidence  Values 


1.  INTRODUCTION 

Observations  of  the  total  electron  content  (TKC)  of  the  ionosphere  have  been 
made  over  the  past  two  decades  and  have  been  used  extensively  to  improve  our 
knowledge  of  the  ionosphere  and  to  correct  for  its  effects  on  signals  propagating 
through  it.  Most  observations  of  TEC  are  made  using  a  satellite  not  directly  over¬ 
head  of  the  observer,  yielding  what  can  be  called  TEC  ,  the  total  electron  content 
along  a  slant  path.  To  bo  more  readilv  understood  and  more  useful  in  studies  of 
the  ionosphere  and  for  updating  models  of  the  ionosphere  for  operational  purposes, 

these  observations  must  be  converted  to  a  vertical  incidence  value,  TEC  ,  corres- 

v 

ponding  to  a  vertical  profile  at  some  point  in  the  ionosphere  along  the  ground- 
sc  ellite  path.  This  conversion  is  bv  necessity  an  approximate  procedure,  the  limits 
of  which  will  be  described  in  this  report. 

An  approximate  relationship  between  TEC's  and  TKC'  mav  be  established  bv  a 
little  geometry,  as  illustrated  in  Figure  1.  We  have 
HS 

TKC  -  /  N  ds  tit 

0 
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Figure  1.  Observer -Satellite  Geometry.  The  satellite  is 
at  an  altitude  h;  the  elevation  angle  at  the  ground  is  Id; 
and  the  angle  of  incidence  at  the  satellite  is  I 


where  the  integration  is  performed  along  the  'observer -satellite  path  and  HS  is  the 
height  of  the  satellite.  Using  Figure  1,  it  follows  that 
HS 

TEC  =  /  Nseeldh  (2) 

s  0 


where  the  integration  is  now  performed  in  the  radial  direction  and  1  is  the  angle  of 
incidence  of  the  ray.  An  expression  for  I  can  be  deduced  from  Figure  1.  using  the 
sine  rule: 


secKh,  E)  - 


R  cos  E 
He  h 


-1/2 


(3) 


where  E  is  the  angle  of  elevation  at  thi  surface  of  the  earth. 

With  the  loss  of  some  accuracy,  sec  I  may  be  removed  from  the  integrand  and 
replaced  in  Eq.  (2)  bv  an  average  value,  seTf; 


(4) 


HS 

TEC  =  iecT  /  N  dh  . 
s  0 

The  integrand  on  the  right  hand  side  of  Eq.  (3)  is  just  the  vertical  incidence  total 

content,  TEC  ,  so  we  have 
v 

TEC  =  secT  TEC  .  (5) 

s  v 

This  is  the  required  equation  relating  slant-  and  vertical-incidence  values  of  TEC. 

The  use  of  an  average  value  of  sec  I,  with  sec  1  being  evaluated  at  an  appro¬ 
priate  altitude  hj,  implicitly  assumes  that  the  integrand  N(h)  sec  I  (h.  E)  peaks  fairly 
sharply  at  hj.  Since  sec  1  is  a  slowly  varying  function  of  h,  the  integrand  in  fact 
peaks  around  the  peak  of  the  .N(h)  uistribution  and  hj  should  be  set  to  an  altitude 
near  hmi'2.  The  error  resulting  from  the  use  of  Eq.  (5)  will  depend  to  an  important 
extent  on  the  choice  of  h.. 

Traditionally,  most  observations  of  TEC  have  been  made  using  the  Faraday 
rotation  technique  and  geostationary  satellites,  l  or  these  observations,  the  ob¬ 
served  values  of  TEC  correspond  to  the  total  content  up  to  2000  km,  and  are  applied 
to  the  ionospheric  point  where  the  ground-satellite  path  cuts  an  altitude  of  400  to 
420  km.  This  procedure  was  recommended  by  Tithe  ridge.  *  The  altitude  of  420  km, 
the  "mean  field  height",  is  appropriate  under  such  conditions  because  the  amount 
of  Faraday  rotation  depends  on  the  earth's  magnetic  field,  as  well  as  on  the  ionos¬ 
phere  and  geometry  of  the  ray  path. 

The  same  choice  of  altitude  is  not  necessarily  appropriate  for  the  conversion 
of  TEC  to  TEC  when  different  measurement  techniques,  not  affected  bv  the 
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earth's  magnetic  field,  are  employed.  Llewellyn  and  Bent,"  for  example,  set  hj 
equal  to  hmF2,  the  height  of  the  1  2  layer  peak,  for  their  model  studies. 

For  reasonably  high  angles  of  elevation  (say  30  to  :(0°),  the  exact  altitude  at 
which  sec  I  is  calculated  (and  which  is  used  to  determine  the  ionospheric  point  with 
which  TEC  is  to  be  associated)  is  not  important  because  TEC  is  not  a  sensitive 
function  of  this  altitude.  However  for  very  low  angles  of  incidence,  such  as  can  be 
encountered  using  polar-orbiting  satellites  and  high- latitude  ground  stations  looking 
at  geostationary  satellites,  sec  I  varies  rapidly  with  height  and  the  value  of  TECv 
determined  from  a  value  of  TEC  depends  strongly  on  the  choice  of  height  at  which 
sec  1  is  evaluated. 


1.  Tithe  ridge,  .1.  E.  (1072)  Determination  of  ionospheric  electron  content  from  the 

Earadav  rotation  of  geostationary  satellite  signals.  Planetary  and  Space  Science, 
20(3):  353-370. 

2.  I.lewellvn,  S.K.,  and  Bent,  H.  B.  i  11*73)  Documentation  and  Description  of  the 

Bent  Ionospheric  Model,  A  Ft '  HI. -Til- 73- 0  K  Is  7 ,  AD  772733. 
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A  second  problem  encountered  with  low  elevation  angles  is  that  the  trans- 
ionospheric  path  is  so  long  that  horizontal  gradients  become  important  and  the 
simplifying  assumptions  on  which  Eq.  t5>  are  based  no  longer  apply. 

This  report  therefore  investigates  the  range  of  validity  of  Eq.  (51  for  low 
elevation  angles.  The  errors  in  Eq.  c 5 *  are  studied  as  a  function  of  station  loca¬ 
tion,  azimuth,  elevation,  local  time,  season,  and  solar  activity,  as  well  as  eleva¬ 
tion  angle.  For  those  regions  in  which  large  errors  are  encountered,  Eq.  (5) 
should  not  be  used.  The  most  appropriate  altitude  at  which  sec  1  should  be  evaluated 
is  also  studied. 

Section  2  describes  the  numerical  methods  used  in  this  analysis.  Section  •> 
discusses  the  errors  associated  with  the  use  of  different  values  of  hj  and  the  most 
appropriate  values  to  choose.  The  conversion  errors  obtained  using  Eq.  (5)  are 
discussed  in  Sections  4.  5.  and  6  for  typical  low,  high,  and  mid- latitude  stations, 
while  Section  7  dismisses  the  errors  for  a  chain  of  stations  at  75'  \V.  The  cases 
studied  are  not  exhaustive,  and  are  meant  to  be  representative  only.  More  detailed 
and  specific  results  may  be  obtained  b\  running  the  compute  -  program  SI  .A  NT  EC 
which  is  available  from  the  author.  The  results  are  summarized  in  Section  8. 


1.  METHOD  OK  ANALYSIS 


In  order  to  detei  mine  the  errors  in  TEC  .  to  TEC.,  conversion,  it  is  necessary 

to  have  accurate  values  lor  each,  in  general,  this  can  be  done  only  by  model  studies, 

since  simultaneous  observations  of  TEC  and  TEC  are  very  difficult  to  arrange. 

1'he  model  of  the  ionosphere  which  is  used  must  be  as  faithful  a  rept  sentation  of 

the  real  ionosphere  as  possible,  especially  in  regard  to  the  variation  of  electron 

2 

density  with  height.  The  limit  model"  is  such  a  model,  having  been  constructed  so 
as  to  reproduce  TEC  observations  over  a  wide  range  of  conditions.  Consequently, 
this  model  has  been  Used  in  the  present  study. 

The  calculation  of  TEC  for  a  given  location  is  readily  achieved  using  the  tech¬ 
niques  used  by  Bent— calculate  foE2  and  hmE2  using  spherical  harmonic  expansions: 
determine  the  corresponding  vertical  Nth)  profile  using  an  empirical  model;  inte¬ 
grate  the  (analytic)  vertical  Nth)  profile  to  obtain  TEC  . 

The  calculation  of  TEC  ,  on  the  other  hand,  is  more  complex  because  the  inte¬ 
gration  is  no  longer  analytic  and  must  be  replaced  by  a  numerical  integration.  We 
have 


US 

TEC  =  /  N(h,  0  )  sec  Ilh,  E>  dh 

s  0 


(6) 
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where  (a)  6  denotes  that  N  can  change  with  latitude  and  longitude,  (b)  I  is  the  angle 
of  the  incidence  at  the  height  h,  corresponding  to  the  elevation  angle  E,  and  (c)  the 
integration  is  performed  from  the  ground  (h=  0)  to  the  satellite  (h=  HS). 

Equation  (61  can  be  replaced  by 

M 

TEC  =  Zl  N.  sec  I.  w.  (7) 

s  ,  1  11 

that  is,  by  an  M-point  quadrature  with  the  integrand  being  evaluated  at  the  heights 
Ir  and  with  quadrature  weights  w. .  The  values  of  h.  depend  on  the  quadrature 
method  and  its  order,  and  also  on  the  limits  of  the  height  integration.  We  have 

3 

chosen  to  use  Gauss- I.egondre  quadrature  (for  example,  Abramowitz  and  Stegun) 
of  different  orders.  The  quadrature  points  and  weights  for  this  method  correspond 
to  integration  limits  of  -1  to  +  1  so  the  method  must  be  generalized  for  use  with 
limits  a  to  b.  This  is  achieved  bv  changing  the  integration  variable: 
b 

[  =  /  fl.vl  dy  .  (8) 


Define 


(9) 


then 


1 


-  1 

/  f |  v  ( x ) )  dx 

-  1 


+  I 

f  f|vtxt]  dx 
-  I 

M 

w.fl.vlxd]  . 


(10) 


The  x.  and  w.  are  the  Gaussian  quadrature  points  and  weights,  y(xd  is  evaluated 
using  Eq.  UD  and  f  is  evaluated  at  yix.l.  In  the  present  example,  f  =  Nsec  I. 


3.  Abramowitz,  IW. ,  and  Stegun,  I.  A.  (1964)  Handbook  of  Mathematical  Functions. 
N11S  Applied  Mathematics  Series  No.  55,  I’.S.  Govt,  r’rinting  Office, 
Washington,  D.  C. 
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In  practice,  the  value  of  TEC  obtained  by  numerical  integration  will  have 
some  error,  the  size  of  which  will  depend  on  the  order  of  integration  chosen.  The 
error  will  decrease  as  the  order  increases,  but  the  computation  time  will  increase, 
so  it  is  necessary  to  choose  as  low  an  order  as  possible,  commensurate  with  the 
desired  accuracy. 

The  required  order  of  integration  may  be  reduced  bv  working  in  terms  of  the 
error  itself.  This  error  is  defined  as 


ERROR  =  seel  •  TEC  -  TEC 

V  V  s 


and  can  be  written  as 


141  141 

ERROR  =  see  1^.  w.(Nv).  -  w. (Nsec  I). 


7,  w.  {seel  '  (N  * 

:  ,  1  I  4  4  1 


•  1 
1  I 


In  the  first  term,  seel  is  a  constant.  X  is  the  vertical  profile  pissing  through 
the  point  at  tvliieh  see  1  is  calculated,  and  is  evaluated  at  ea<  h  of  the  M  quadrature 
heights.  In  the  second  term,  X  is  the  vertical  profile  passing  through  each  quadra¬ 
ture  point  along  the  integration  path,  seel  and  X  are  both  evaluated  at  each 
quadrature  height. 

The  required  order  of  integration  may  be  further  reduced  bv  noting  that  the 
error  term  changes  sign  at  the  height  at  which  seel  is  evaluated.  Consequently, 
it  is  more  accurate  to  use  two  ranges  of  integration.  111  to  hj  and  h.  to  H\  ,  where 
III.  and  IIL  are  tin-  lower  and  upper  limits  of  the  full  integration. 

Once  the  error  term  has  been  calculated,  the  correct  value  of  TEC  can  be 

s 

obtained  from  Eq.  (11).  Most  of  the  discussion  of  errors  in  this  report  is  con¬ 
cerned  4Vith  the  error 


ERROR 

TEC 


This  error  is  positive  if  the  simple  see  I  approximation  yields  too  high  a  result, 
and  vice  versa. 

1'igure  2  shows  some  typical  integrands  (terms)  of  Eq.  (12).  It  can  be  seen 
that  the  error  changes  sign  at  hj,  which  is  420  km  in  this  ease.  The  limits  of 


integration  are  200  and  2000  km  and  8-point  quadrature  is  used  in  both  integration 
ranges.  Since  the  error  is  an  oscillating  function  of  height,  some  cancellation 
occurs,  the  over-all  error  depending  on  the  extent  of  the  cancellation.  The  upper 
limit  of  integration  could  be  lowered  to  1200  km  without  introducing  any  appreciable 
errors,  since  the  contribution  to  the  integral  from  1200  to  2000  km  is  negligible. 


Figure  2.  Variation  of  the  Error  Integrand  With  Integration 
Height  for  Three  Local  Times  (Kwajalein,  March,  11=  100, 
hj  =  420  km.  Elevation  =  10°) 
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3.  SELECTION  OF  THE  HEIGHT  h( 

The  errors  in  the  conversion  of  TEC  to  TEC  are  found  to  be  a  complicated 

s  v 

function  of  ali  variables  — station  location,  solar  activity,  month,  time  of  day, 
azimuth,  elevation,  hj— and  these  variations  are  inter-dependent.  Consequently, 
in  order  to  investigate  each  different  variation  it  is  necessary  to  hold  all  other 
parameters  constant  or  to  average  the  errors  over  the  ranges  of  these  other 
parameters. 

The  first  variation  which  must  be  considered  is  that  of  the  errors  versus  the 
choice  of  hj.  Ideally,  a  single  value  of  hj  should  lead  to  the  smallest  possible 
errors  under  all  circumstances.  However,  this  is  found  to  be  not  the  case. 

The  conversion  error  has  been  studied  as  a  function  of  hj,  for  three  stations, 
four  months  of  the  year,  and  for  low  and  high  solar  activity.  The  stations  chosen 
are  representative  low-latitude  (Kwajalein),  mid-latitude  (Cape  Canaveral),  and 
high-latitude  (Goose  Bay)  stations.  Values  of  hj  have  been  taken  from  320  to 
4  60  km,  in  steps  of  20  km,  and  values  related  to  hml'2  have  also  been  used;  hmF2 
plus  0,  50.  100  and  150  km.  The  use  of  a  fixed  altitude,  rather  than  one  related  to 
hml'2,  seems  preferable  because  it  obviates  the  need  for  a  knowledge  (most  likely 
from  an  ionospheric  model)  of  hml'2.  On  the  other  hand,  however,  the  appropriate 
value  of  h|  would  be  expected  to  be  related  to  hml'2  and  have  similar  diurnal  and 
other  variations. 

Table  1  shows  the  KMS  errors  for  Kwajalein,  Cape  Canaveral  and  Goose  Bay, 
for  March,  June,  September,  and  December,  and  for  =  100  and  10  (corres¬ 
ponding  fluxes  of  150  and  75).  Only  those  errors  necessary  to  define  the  pattern 
have  been  calculated.  The  errors  have  been  averaged  over  the  day  (00,  06,  12, 

1 H  i; T) ,  azimuth  ( 0" ,  00°,  160°,  270  ’>  and  elevation  angle  (0%  10°,  20°,  30°).  The 
entries  underlined  in  the  tabic  are  the  minimum  values  for  each  case.  Clearly, 
there  is  no  single  height,  h(,  which  should  be  used  in  all  circumstances.  For 
Cape  Canaveral,  the  error  minimizes  at  360  ±  20  km,  while  for  Kwajalein  it 
minimizes  at  400  ±  20  km  (approximately).  For  Goose  Bay,  low  solar  activity, 
the  error  minimizes  at  360  ±  20  km,  but  for  high  solar  activity,  there  is  no  obvious 
best  choice  of  hj.  (This  indicates  the  presence  of  a  systematic  error  in  the  conver¬ 
sion  which  cannot  be  overcome  by  using  a  suitable  hj>. 

When  hj  is  related  to  hmF2,  the  errors  are  least  for  hmK2  +  50  or 
hml'2  +  100  km,  but  have  not  been  investigated  in  great  detail.  This  result  is 
readily  explained  by  noting  that  the  centroid  of  the  Bent  X(h>  profile  is  about  50  km 
above  hml'2.  Two  things  can  be  noted  about  the  results  for  h.  related  to  hmF2: 


Table  1.  RMS  Percentage  Conversion  Errors  Averaged  Over  Elevation  Angle 
(0,  10,  20,  30°),  Azimuth  (0,  90,  180,  27  0° )  and  Over  the  Day  (0,  06.  12, 

18  UT).  The  errors  are  given  for  different  values  of  hj,  four  months,  two 
levels  of  solar  activity  and  for  stations  at  three  latitudes.  The  underlined 
values  are  the  minimum  errors  as  a  function  of  hj,  for  each  set  of  h, 

(hm  F2  •  0.  50,  100,  150  km;  320 -460  km )  1 


Station 

Month 

R 

HM 

50 

100 

150 

320 

340 

360 

380 

400 

420 

440 

460 

Kwajalein 

MAR 

100 

9.  3 

6.  3 

6.  3 

8.  0 

12.  1 

10.  4 

9.  2 

8.  4 

7.9 

7.7 

7.  8 

JIN 

9.5 

5.  8 

5.  3 

7.2 

12.  7 

10.  9 

9.  3 

8.  2 

7.5 

7.  1 

7.  1 

sep 

8.  8 

5.4 

5.  1 

7.  1 

10.  5 

9.  1 

8.0 

7.2 

6.  8 

6.  9 

7.  2 

DEC 

8.  8 

4.  8 

5.  8 

8.  6 

10.  9 

9.  0 

7.  6 

6.  8 

6.  8 

7.  1 

7.7 

MAR 

10 

10.  7 

6.  4 

5.  8 

8.  0 

7.  5 

7.0 

6.  9 

7.  1 

7.  8 

8.  6 

9.  6 

JIN 

1  1.  9 

7.0 

4.  9 

6.  1 

7.7 

6.  5 

5.  9 

5.  9 

6.  3 

6.  9 

7.8 

SEP 

12.  1 

7.0 

5.  6 

7.  8 

8.  0 

7.  0 

6.  6 

6.  8 

7.  6 

8.  5 

9.  6 

DEC 

10.  3 

5.  4 

4.  4 

7.  0 

6.  5 

5.  4 

5.  0 

5.  3 

6.  1 

7.  2 

8.  2 

Cape 

MAR 

100 

6.  0 

2.  8 

3.  6 

6.  0 

5.  6 

4.  2 

3.  2 

2.  9 

3.  2 

4.  0 

4.  9 

Canaveral 

J  IN 

6.  5 

3.  3 

3.  3 

5.  5 

7.3 

5.  7 

4.  2 

3.  0 

2.  3 

2.  4 

3.  1 

SEP 

6.  0 

4.  5 

3.  3 

2.7 

2.8 

3.  5 

4.4 

DEC 

3.  6 

2.  9 

3.  0 

3.  9 

4.  9 

6.  1 

MAR 

10 

9.  n 

4.8 

2.  1 

4.  1 

5.7 

4.  1 

2.  9 

2.  5 

2.  9 

3.  9 

4.  9 

JIN 

9.  4 

5.  1 

2.  7 

3.  7 

6.  1 

4.  6 

3.  4 

2.  6 

2.  5 

3.  0 

3.  8 

SEP 

6.  2 

4.  6 

3.  3 

2.7 

2.  9 

3.  7 

4.  7 

DEC 

3.  6 

2.  7 

2.  8 

3.  6 

4.  7 

5.  8 

Goose  Bay 

MAR 

100 

12.  1 

3.  4 

4.  0 

3.  5 

11.  5 

9.  5 

7.8 

6.  2 

5.  0 

4.  2 

3.  8 

JEN 

10.  0 

8.  4 

6.  9 

5.  5 

4.  3 

3.  3 

2.5 

SEP 

9.  4 

7.  7 

6.  2 

4.  9 

3.  8 

3.  2 

3.  0 

DEC 

12.  6 

10.  2 

8.  2 

6.  6 

5.  4 

5.  0 

5.  2 

MAR 

10 

8.  7 

4. 2 

2.  2 

4.  1 

3.  6 

2.  4 

1.  9 

2.4 

3.  3 

4.  3 

5.  3 

JEN 

4.  3 

3.  0 

1.  9 

1.  6 

2.  1 

3.  0 

4.  0 

SEP 

3.  2 

2.  2 

2.  0 

2.  5 

3.  4 

4.4 

5.4 

DEC 


2.4  2.2  2.9  4.0  5.  1  6.2  7.3 


(1)  For  Kwajaiein,  at  least,  the  RMS  errors  are  lower  for  hj  =  hmF2  +  100  km 
than  for  hj  =  400  km.  This  indicates  that  the  former  would  be  preferable,  provided 
that  hmF2  were  known  reliably. 

(2)  The  choice  of  hmF2  +  50  or  +  100  km  is  usually  far  superior  to  the  use 

2 

of  hmF2  itself,  which  is  the  procedure  adopted  by  Llewellyn  and  Bent. 

The  diurnal  variation  of  the  error  for  values  of  hj  of  400  km  and  hmF2  +  50  km 
is  shown  in  Figure  3,  for  four  azimuths.  The  results  are  for  Kwajaiein.  March, 
l(  --  100,  elevation  angle  equal  to  10°.  The  largest  errors  are  obtained  looking 
north  and  south  (the  reasons  for  this  are  explained  in  Section  4)  and  are  greater  in 
general  for  hj  =  400  km.  These  errors  tend  to  be  positive  during  the  day  (hj  is  too 
low)  and  negative  at  night  <hj  too  high).  When  hj  is  allowed  to  follow  the  variation 
of  hmF2,  the  errors  generally  decrease  below  those  for  h^  =  400  km  and  to  have  a 
less  definite  diurnal  variation. 


Figure  3.  Diurnal  Variation  of  the  Conversion  Error  for  North,  South,  East, 
and  West  Aziniuths,  for  h.  =  400  km  and  hmF2  +  50  km  (Kwajaiein,  March, 
R=  100,  Elevation  =  10") 
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Figure  4  illustrates  how  the  error  depends  on  the  choice  of  hj,  for  four  azi¬ 
muths  and  four  local  times.  The  error  is  plotted  as  a  function  of  elevation  angle 
and  the  results  are  for  Kwajalein.  March,  R  =  100.  It  can  be  seen  that  the  error 
is  a  complicated  function  of  azimuth,  elevation,  local  time,  and  hj,  at  least  for  low 
elevation  angles,  but  it  is  the  variation  of  error  with  hj  which  is  of  immediate  con¬ 
cern  here.  This  is  especially  important  for  low  elevation  angles,  whereas  for 
higher  elevation  angles  (greater  than  20°  or  30°)  the  errors  are  fairly  insensitive 
to  the  choice  of  hj.  Only  for  low  elevation  angles  is  it  really  necessary  to  attempt 
to  minimize  the  errors  by  choosing  the  "best"  value  of  hj. 
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Figure  4.  Variation  of  the  Conversion  Error  With  Elevation  Angle  for  North, 
East,  South,  and  West  Azimuths,  for  Four  Local  Times  (1  1,  17,  23,  05),  and 
lor  Three  Values  of  h.,  (360,  400  and  440  km)  (Kwajalein,  March,  R  100) 


The  other  variations  implicit  in  Figure  4  are  discussed  in  the  next  section. 

The  results  discussed  in  the  remainder  of  this  repori  were  obtained  using 
hj  =  400  km.  This  seems  a  good  average  of  the  best  values  and  is  an  appropriate 
height  for  F- region  studies.  For  the  results  presented  in  Table  1.  the  average 
error  is  4.7%  for  hj  =  400  km.  It  should  be  noted  that  the  errors  do  not  increase 
dramatically  for  other  values  of  hj  —  for  the  results  presented  in  Table  1,  the  errors 
increase  only  to  5.  5%  for  hj  =  -160  and  440  km. 

4.  CONVERSION  ERRORS  FOK  \  I  OW-LVTITl  l)E  ST  \TION 

Kwajalein  (4.  3°  N,  167.  r  E,  8.  6^  dip!  lies  between  the  crests  of  tile  equatorial 
anomaly  and  can  be  regarded  hereasa  typical  low-latitude  station.  The  variation 
of  the  conversion  errors  for  March,  R  -  100  have  been  given  in  Figure  4,  for  four 
azimuths  and  four  local  times.  The  errors  are  clearly  greatest  when  looking 
north  and  south  at  low  elevation  angles,  f  igure  a  shows  the  errors  for  March, 

June,  and  December,  for  hj  -  400  km,  and  for  R  100,  four  azimuths  and  four 
local  times.  The  errors  are  again  greatest  when  looking  north  .  nu  south,  the 
greatest  values  being  for  looking  south  in  March  during  the  day  ,  e-  ,  teding  20'  . 

As  might  be  expected,  the  errors  are  largest  when  there  >re  1  .  rg  r  horizontal 
gradients  of  electron  density  along  the  observer-satellite  path.  I  'gure  t>  -  hows 
the  horizontal  ranges  (in  degrees)  at  which  a  ray  from  the  satellite  will  cut  the 
Ionosphere  at  a  specified  altitude,  for  elevation  angles  of  (l  to  30  .  Propagation 
is  either  north-south  or  south-north  for  thus  example  since  the—  azimuths  corres¬ 
pond  to  the  greatest  errors.  A  ray  at  zero  degrees,  tor  example,  will  travel  about 
20’  in  latitude  before  encountering  the  iono.-phero  at  400  kin.  I  iguros  7a  and  7b 
show  how  the  ionosphere  varies  along  the  integration  path  for  March,  R  -  i  00, 
hj  --  400  km,  elevation  =  10°,  local  time  -  1100,  which  is  when  flu  i  (inversion  error 
has  some  of  its  highest  values.  It  can  he  seen  that  under  these  conditions  the  ray 
passes  right  through  a  large  part  of  the  crests  of  the  equatorial  anomaly,  indicating 
that  it  is  the  presence  of  the  anomaly  which  cause  the  large  errors.  The  errors 
are  significantly  lower  for  the  east  and  west  directions  (Fig tin-  ">l. 

Evidence  supporting  the  contention  that  it  is  horizontal  gradients  which  cause 
the  large  errors  is  shown  in  Figure  H.  which  shows  plots  of  percentage  error  versus 
gradient  for  the  four  azimuths  (March,  R  =  100,  elevation  10  ,  hj  400  km'.  The 
gradients  are  represented  by  the  change  of  liml-2  ..nd  fol  2  over  the  height  range  of 
the  first  integral  (200  to  400  km).  Each  data  point  is  for  a  particular  time  of  day 
from  0  to  22  I  T,  in  steps  of  two  hours.  The  two  boxed  values  correspond  to  the  data 
plotted  in  f  igures  7a  and  7b,  and  the  straight  lines  are  the  least  squares  fit  lines  of 
error  on  gradient.  The  correlation  coefficients  are  very  high  for  north  (0.  07)  and 
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Figup  H.  Variations  of  the  Kwajalein  Conversion  Errors  With  Gradients  in 
h:nF2  i  top  panel)  and  fol'2  (bottom  panel),  for  North,  South,  East,  and  West 
Azimuths.  The  straight  lines  are  the  least  squares  fit  regression  lines  and 
correspond  to  the  linear  correlation  coefficients,  r.  (March,  H=  100, 
Elevation  •  10  t 


I  he  graphs  >r  Figure  H  support  the  contention  that  the  conversion  errors  are 
due  :•>  gradients  along  the  path  if  hm  F2  and  foF2,  especially  the  former.  The 

■  orrelati  >n  is  most  pronounced  for  north  and  south  azimuths.  The  correlation  with 
hr -i  I'2  rmv  he  explained  qualitatively  bv  noting  that  if  hmE2  increases  as  the  integra- 
'i  in  lieight  increases,  the  same  Nth)  profile  will,  for  a  given  altitude  below  hmF2, 

■  icld  lower  densities  fur  die  higher  values  of  hmF2.  It  is  difficult  to  be  specific 
about  die  exact  source  of  the  conversion  errors,  because  the  net  error  is  a  result 
of  cancellation  between  the  errors  above  and  below  hj,  and  each  of  these  depends 
in  the  gradients  in  foF2  and  hm  F2. 


23 


5.  CONVERSION  ERRORS  FOR  A  HIGH-LATITUDE  STATION 


Having  found  that  the  conversion  errors  at  a  low -latitude  station  are  mainly 
due  to  the  presence  of  horizontal  gradients  associated  with  the  equatorial  anomaly, 
it  is  logical  to  choose  as  a  representative  high -latitude  station  one  which  can  "see" 
the  main  feature  of  the  high-latitude  ionosphere,  the  high -latitude  trough.  Conse¬ 
quently  Goose  Bay  (53.  3°  N,  299.  5“  E,  65“  N  geomagnetic)  has  been  chosen  as  a 
representative  high -latitude  station.  Its  position  relative  to  the  trough  under  cer¬ 
tain  conditions  is  shown  in  Figure  9. 


90* 


Figure  9.  Location  of  the  High  Latitude  Trough 
for  Magnetically  Quiet  Conditions  (March,  K  50, 
06  L'T),  After  Halcrow  and  Nisbet^ 


4.  Halcrow,  B.  M.  ,  and  Nisbet,  .I.S.  (1977)  A  model  of  the  l-'2  peak  electron  den¬ 
sities  in  the  main  trough  region  of  the  ionosphere.  Radio  Science, 

1 2  ( 5 ) :  8  15 -820. 

MM 


24 


The  range  of  latitudes  which  can  be  seen  from  Goose  Bay  is  illustrated  in 
Figure  10.  (It  should  be  noted  here  that  the  Bent  model  of  the  ionosphere  is  uncer¬ 
tain  for  very  high  latitudes.  )  Figures  11a  and  11b  show  how  hmF2  and  foF2  vary 
along  the  integration  paths  for  four  azimuths  (March.  R  -  100,  hj  =  400  km,  0°  eleva 
tion,  06  UT).  The  associated  errors  exceed  10%  in  all  cases. 


Figure  10.  Latitudinal  Coverage  of  the  Ionosphere  From 
Goose  Bav,  for  Different  Flevation  Angles  and  Ionospheric 
Altitudes 

The  errors  are  shown  as  a  function  of  azimuth,  elevation  angle,  local  lino 
and  season  in  Figure  12  ( h j  -  400  km).  In  general,  the  major  errors  are  for  cast, 
west  and  south  propagation  at  around  02  1.1  ,  when  the  line  of  sight  cuts  son  e  part 
of  the  trough.  For  north  azimuth,  the  ionospheric  part  of  the  rav  path  lies  north 
of  the  trough  and  the  errors  are  somewhat  smaller.  There  iocs  not  sennit  to  he 
anv  sin  pie  relationship  between  the  errors  and  gradients  in  either  hml  2  >r  tol-2. 
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Figure  12.  Variations  of  the  Goose  Bay  Conversion  Errors  With  Elevation 
Angle  for  North,  East,  South,  anti  West  Azimuths,  for  Four  Local  Times 
(20,  02,  08,  14)  and  for  Three  Seasons  (March,  June,  December), 

<11  100.  h,  =  400  km  i 

(»,  CONVERSION  ERRORS  FOR  A  MID-EATITl  DE  STATION 

Cape  Canaveral  (28.43  N.  279.  4°  E,  39.  6’ N  geomagnetic)  has  been  chosen  as 
a  tvpiral  m id -latitude  station.  The  part  of  the  ionosphere  covered  bv  looking  north 
and  south  from  Cape  Canaveral  is  illustrated  in  Figure  13.  Figures  14a  and  14b 
correspond  to  the  cases  with  largest  errors  (for  March,  R  =  100,  hj  =  400  km). 

The  large  error  obtained  looking  south  is  a  result  of  the  ray  passing  through  the 
northern  crest  of  the  equatorial  anomaly,  while  the  large  error  looking  north  is  a 
result  of  the  lowering  of  the  ionosphere  with  increasing  latitude. 

The  general  correlation  of  the  conversion  error  with  gradients  in  hmF2  and 
loF'2  is  illustrated  in  Figure  15  (March,  R  -  100,  hj  -  400  km,  0°  elevation)  for  four 
azimuths.  The  gradients  are  represented  bv  the  changes  of  hmF2  and  foF2  over  the 
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HEIGHT  ALONG  PATH  (km) 

Figure  13.  Latitudinal  Coverage  of  the  Ionosphere  From 
Cape  Canaveral  for  Different  Elevation  Angles  and  Ionospheric 
Altitudes 

first  integration  range,  while  the  azimuths  are  represented  by  the  letters  N,  E,  S 
and  W.  There  are  four  data  points  (00,  06,  12,  18  UT)  for  each  azimuth.  The  re¬ 
gression  lines  plotted  are  for  all  16  data  points  and  show  a  much  higher  correlation 
of  error  with  the  gradient  in  hmF2  (r  -  0.68)  than  with  foF2  (r  -  0.  15).  Note  that 

hm  F2  changes  bv  onlv  13  km  at  most  (when  looking  south)  in  the  5  degrees  of  latitude 
encompassed  bv  the  first  integration  range.  The  correlation  with  the  gradient  in 
hmF2  is  highest  when  looking  south  (r  --  0.  80)  and  lowest  when  looking  west  (r  =  0.  20). 
The  boxed  data  points  for  south  and  north  azimuth  correspond  to  the  cases  shown 
in  Figures  Ha  and  14b. 

The  conversion  errors  for  Cape  Canaveral  are  summarized  in  Table  2,  which 
shows  the  smallest  elevation  angle  for  which  the  errors  are  less  than  5%  at  00.  06. 

12  and  18  FT  (I..T  FT  •  18.5),  and  the  error  corresponding  to  an  elevation  angle 
of  20’.  Both  of  these  levels  are  arbitrary,  but  reasonable  choices.  The  parameter 
h|  has  been  set  to  400  km,  and  the  results  are  four  azimuths,  four  times  of  day, 
three  seasons,  and  two  levels  of  solar  activity. 
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Figure  1-1.  Variations  oT  foF2  anil  hm K2  Along  ihe  Integration  Path 
Prom  rape  Canaveral  for  l  wo  Times  When  the  Conversion  Krror 
Was  Fairlv  High.  (March.  It  =  100,  h.  400  km) 


figure  15.  Variations  of  the  Cape  Canaveral  Conversion  i’.rrors  With  Gradients 
in  hntl'2  and  foFU,  for  North,  South,  Hast,  and  West  Azimuths.  The  straight 
lines  are  the  least  squares  fit  regression  lines  and  correspond  to  the  linear 
correlation  coefficients,  r.  (March.  H  100,  hr  ~  400  km,  elevation  O’). 

Ihe  boxed  values  correspond  to  the  conditions  of  Figure  14 


Table  2.  Conversion  Errors  for  Cape  Canaveral  Showing  the  Smallest  Elevation 

Angle  for  Which  the  Error  is  Less  Than  5%  at  0,  06,  12  and  18  UT 

(LT  =  UT  +  18.5),  and  the  Error  Corresponding  to  an  Elevation  Angle  of  20° 


Elevation  Angle  at  Which 

Error 

S  5% 

Error  at  20 

0  Elevation 

0 

90 

180 

270 

0 

90 

180 

27  0 

March 

0 

20 

20 

0 

0 

-3 

-4 

-  1 

-1 

R  =  0 

6 

0 

0 

0 

0 

3 

2 

5 

0 

12 

10 

0 

10 

20 

-1 

2 

0 

4 

18 

0 

0 

0 

0 

-2 

0 

2 

0 

R  100 

0 

30 

20 

10 

0 

-5 

-4 

-i 

0 

6 

0 

0 

40 

0 

-1 

2 

6 

1 

12 

0 

0 

10 

0 

-3 

3 

0 

-3 

18 

0 

0 

10 

0 

-2 

1 

-1 

2 

Aver. 

RMS 

2.  8 

2.  6 

2.  9 

2.  0 

2.  6 

June 

0 

0 

0 

0 

0 

-2 

- 1 

2 

- 1 

R  -  10 

6 

0 

0 

40 

0 

5 

0 

7 

1 

12 

0 

0 

0 

0 

1 

1 

3 

0 

18 

0 

0 

0 

0 

1 

1 

0 

2 

R  =  100 

0 

0 

0 

10 

0 

-2 

0 

1 

0 

6 

0 

0 

40 

0 

3 

1 

6 

2 

12 

0 

0 

40 

0 

- 1 

1 

5 

i 

18 

0 

0 

0 

0 

- 1 

2 

i 

2 

Aver. 

RMS 

2.  4 

1. 1 

4.  0 

1.  4 

2  2 

Dec 

0 

10 

20 

0 

0 

0 

-4 

0 

2 

R  =  10 

6 

0 

0 

10 

0 

2 

2 

3 

0 

12 

20 

0 

10 

30 

-3 

2 

1 

-5 

18 

10 

20 

20 

10 

-2 

-3 

-3 

-2 

R  100 

0 

20 

20 

0 

0 

-4 

-4 

_  9 

i 

6 

0 

0 

0 

0 

-4 

2 

4 

i 

12 

20 

30 

10 

20 

-4 

5 

2 

-4 

18 

20 

10 

30 

10 

-4 

-2 

-6 

-  1 

Aver. 

RMS 

3.  2 

3.  2 

3.  1 

2.  5 

3.0 

Aver, 
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Consider  first  the  elevation  angles.  A  count  of  the  number  of  times  each 
elevation  angle  occurs  in  the  tabLe  shows  that  elevation  angles  of  40°  are  required 
onlv  when  looking  south  (azimuth  =  180°),  and  that  all  four  cases  occur  at  06  or 
12  UT  (0030,  0630  LT).  In  the  opposite  sense,  a  zero  elevation  angle  produces 
errors  less  than  5%  on  most  occasions  when  looking  west. 

The  variation  of  error  with  azimuth  is  better  studied  using  the  errors  corres¬ 
ponding  to  an  elevation  angle  of  20°,  also  shown  in  Table  2.  RMS  errors  have  been 
calculated  for  each  azimuth  and  month  (N  =  8)  and  then  averaged  across  the  azimuth 
values  for  each  month  and  across  the  months  for  each  azimuth.  It  can  be  seen  that 
the  errors  are  greatest  when  looking  south;  greatest  during  December  (winter); 
least  when  looking  west;  and  least  during  June  (summer).  In  general,  the  errors 
are  lower  when  looking  east  or  west  than  when  looking  north  or  south. 

The  conversion  errors  are  also  illustrated  in  Figure  16,  which  shows  the 
variation  of  error  with  elevation  angle  for  north,  east,  south  and  west  azimuths, 
for  four  local  times  (1!),  01,  07,  13),  and  for  three  seasons  (March,  June  and 
December).  The  errors  do  not  show  the  larger  values  attained  near  the  equatorial 
anomaly  (Figure  5)  and  high  latitude  trough  (Figure  12)  (R  100,  hj  =  400  km). 


7.  (X)NV  KKSION  ERRORS  FOR  STATIONS  AT  75°  W 

This  section  describes  in  a  broad  fashion  the  variation  of  the  conversion  error 
for  stations  at  283°  K  (75’  W )  and  at  latitudes  of  -10“  to  70°  (geomagnetic  latitudes 
of  '1°  and  '  82’ N).  The  results  thus  include  the  effects  of  the  high -latitude  trough 
and  the  equatorial  anomaly,  as  well  as  illustrating  the  errors  encountered  when 
neither  of  these  ionospheric  features  need  be  considered.  The  results  are  all  for 
h j  400  km  and  R  100,  and  are  given  separately  for  four  azimuths.  Table  3,  for 
example,  gives  the  results  for  March,  and  includes  (a)  the  elevation  angle  required 
to  be  exceeded  in  order  for  the  RMS  error  (averaged  over  the  eight  times 
01,  04  .  .  .  22  LT)  to  be  %  3%,  or  for  all  eight  errors  to  be  %  5%  in  absolute  values. 
Elevation  angles  between  20’  and  50’  inclusive  are  considered;  (b)  the  RMS  errors 
at  20°  elevation.  This  elevation  is  a  likely  operational  minimum  set  by  antenna  con¬ 
siderations;  (c)  the  RMS  errors  at  30°  elevation.  This  elevation  would  allow  stations 
10*  apart  to  have  overlapping  fields  of  view. 

Tables  4  and  5  give  the  corresponding  results  for  June  (summer)  and  December 
(winter).  The  three  tables  present  no  surprises  and  the  broad  picture  is  the  same 
as  that  deduced  from  the  results  for  Kwajalein,  Cape  Canaveral  and  Goose  Bay: 

(a)  the  errors  are  least  when  looking  east  or  west,  except  at  high  latitudes,  where 
the  trough  is  encountered;  (b)  for  north  and  south  azimuths,  the  errors  are  greatest 
when  the  trough  and  equatorial  anomaly  are  encountered;  (c)  the  general  level  of 
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Figure  16.  V  ariations  of  the  Cape  Canaveral  Conversion  Errors  With  Elevation 
Angle  lor  North.  East.  South,  and  West  Azimuths,  for  Four  Loral  Times 
(10.  01,  07,  El)  and  for  Three  Seasons  (March,  dune,  December),  (R  =  100, 
lij  -  400  km ) 


error  is  greatest  in  December  (winter)  and  least  in  dune.  Quantitative  estimates  of 
the  errors  mav  be  obtained  from  the  tables. 

The  use  of  RMS  errors  averaged  over  the  eight  local  times  hides  any  diurnal 
variation  of  the  error  which  at  times  is  very  significant.  Accordingly,  the  diurnal 
variation  of  the  error  has  been  plotted  in  Figures  17,  18  and  19  for  each  azimuth 
and  station  latitude  and  for  March,  dune  and  December.  The  plots  also  show  up 
anv  tendency  for  a  systematic  error  to  occur.  The  results  shown  are  for  R  =  100, 
hj  400  km  and  an  elevation  of  20’. 


i 
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unction  of  Latitude  and  Azimuth,  (R  -  100,  h.  400  km.  Elevation 


The  most  significant  diurnal  variations  occur  for  low -latitude  observations 
affected  by  the  equatorial  anomaly  (looking  north  or  south  in  any  month).  Significant, 
but  somewhat  smaller  variations  are  found  for  nigh-latitude  stations  affected  by  the 
high -latitude  trough  (for  example,  looking  east  or  west  in  December).  Systematic 
errors  occur  for  a  station  at  10°  S  (geomagnetic  latitude  ~0°)  when  looking  south 
in  March  and  June,  and  in  general  for  stations  north  of  about  50°.  Some  of  the 
systematic  errors  could  be  removed  bv  a  different  choice  of  hj. 

A  third  source  of  gradients  which  might  be  expected  to  yield  significant  errors 
are  the  sunrise  and  sunset  transitions,  especially  the  former.  Such  errors  are 
not  obvious  in  the  results  although  there  are  suggestions  of  dawn  efiects  on  some 
diurnal  plots.  For  example.  Figure  3  for  hj  =  400  km  does  show  a  local  change  of 
error  between  0100  and  0600  l.T  when  looking  east.  Similar  local  increases  can  be 
seen  in  Figures  17  (March)  and  10  (December)  for  station  latitudes  of  -10°  to  30°, 
looking  cast  at  0500  FT,  hut  not  in  Figure  13  (June). 

The  rapid  morning  increase  of  foF2  leads  to  significant  gradients  on  paths 
looking  west,  resulting  at  times  in  noticeable  errors.  For  example,  the  calculated 
values  of  TFC  based  on  die  vertical  incidence  value  are  too  low  near  08  LT  in  March 
and  December  (Figures  17  and  1!*).  live  errors  result  from  a  decrease  in  foF2  and 
an  increase  in  hml-2  along  the  integration  path,  which  looks  towards  earlier  time 
/ones. 


ft.  DISCI  SSION 

The  errors  m  the  conversion  of  slant  TFC  to  equivalent  vertical  incidence 
values  have  been  found  to  varv  with  the  elevation  angle  to  the  satellite,  with  the 
altitude  at  which  seel  is  evaluated,  with  the  presence  of  ionospheric  gradients  along 
tile  ground -satellite  path,  and  with  the  well-known  variations  of  the  ionosphere 
itself  — diurnal,  seasonal,  solar  evade,  and  location.  If  the  elevation  angle  is  re¬ 
stricted  to  20°  or  greater,  the  RMS  errors  averaged  aver  a  day  are  usually  less  than 
about  5%  except  where  the  rav  path  cuts  either  the  equatorial  anomaly  or  the  high- 
latitude  ionospheric  trough.  In  general,  the  largest  errors  are  encountered  near 
the  equatorial  anomalv,  where  they  can  exceed  15%  for  an  elevation  angle  of  20°. 

If  elevation  angles  lower  than  20  ’  are  considered,  conversion  errors  in  excess.of 
25%  are  possible  near  the  anomalv,  while  errors  of  up  to  20%  are  possible  near  the 
high -latitude  trough. 

There  are  several  practical  considerations  associated  with  the  errors  in  the 
c  onveision  process.  Consider  firstly  the  situation  of  mid-latitude  stations  moni¬ 
toring  beacons  on  geostationary  satellites  at  other  longitudes.  In  such  cases,  the 
elevation  angle  can  get  verv  low.  possibly  resulting  in  non -negiligible  conversion 
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errors.  If  the  elevation  angle  gets  down  to  10°,  a  typical  conversion  error  is  5 
(Cape  Canaveral,  Figure  16),  but  can  reach  10%  for  a  fairlv  high  mid-latitude  station 
(Goose  Bay,  Figure  12,  south  azimuth).  The  error  has  a  diurnal  variation,  which 
must  also  be  taken  into  account.  Changes  of  TFC  less  than  5  to  10%  are  not  sig¬ 
nificant  in  such  cases. 

The  next  point  to  consider  is  the  absolute  size  of  the  errors  at  different  loca¬ 
tions  when  the  relative  errors  are  the  same.  A  20%,  error  at  Goose  Hay  at  02  FT 
in  December  (  K  =  100)  is  equivalent  to  a  change  of  TFC  of 1 .  :(  ■  lO16  m  ", 
whereas  the  same  percentage  error  at  Kwajalein  at  17  I  T  in  March  ill  -  100)  is 
equivalent  to  a  much  larger  change  of  ~  8  >  101U  m  “F  In  general,  the  absolute 
errors  at  low-latitude  stations  are  higher  than  at  high-latitude  stations  because  of 
the  higher  electron  densities  involved. 

The  errors  derived  for  low-  and  high-latitude  stations  apph  onh  to  tile  aver¬ 
age  behavior  of  the  ionosphere  under  magnetically  quiet  conditions.  1  ruler  disturbed 
geomagnetic  conditions,  the  high-latitude  trough  expands  equatorwa  rds  and  a  mid- 
latitude”  station  can  come  under  the  trough  and  become  ,  high-lutitud.  station  as 
far  as  TfciC  conversion  errors  are  concerned.  Similarly,  under-  disturbed  geomag¬ 
netic  conditions,  the  poleward  expansion  of  the  crests  of  the  equatorial  :,:iom  ,p. 
is  inhibited  and  a  low-latitude  TFf  station  will  sample  a  sigmiic  nith  different 
section  of  the  anomaly.  Since  it  is  not  possible  to  model  ttie.se  rid  si:  . , !  1  >  r  .j  -to¬ 
day  variations  of  the  ionosphere  and  its  gradients,  the  conversion  error.-  i 

with  a  location  likely  to  become  affected  bv  such  v.,ri  ,tr»n-  -noulti  hi  vioscn 
pessimistically  from  those  given  in  this  report.  In  effect,  tne  stations  should  he 
"moved"  by  cn  amount  comparable  to  the  amount  by  whicti  the  !e  .tares  ot  ttie  ionos¬ 
phere  move.  For  the  equatorial  regions,  this  is  about  t  n  .  I  or  nigh- latitude 
stations,  account  should  be  taken  of  how  far  equatorwa  rd>  the  hign-latr.u  ie  trougli 
has  expanded,  which  will  depend  directly  on  the  level  ot  magnetic  utr.  tty.  4  The 
trough  may  in  fact  be  recognizable  on  the  Flit  records  if  tae.-o  vield  sufficient 
spatial  coverag  . 

The  results  presented  in  this  report  or  obtained  using  further  runs  of  the  pro¬ 
gram  SI.AN'TFC  may  be  used  to  estimate  the  ground-station  separation  required  to 
yield  acceptably  small  conversion  errors.  This  spacing  will  change  with  latitude, 
and  will  also  depend  on  what  level  of  error  is  •.  ceptable.  A  point  which  should  he 
remembered  is  that  even  if  the  conversion  errors  are  relatively  large,  say  20’:.,  it 
may  be  better  to  make  use  of  the  data  rather  than  simply  rejecting  them,  because 
they  do,  after  all,  contain  some  useful  inform;  tion.  In  many  cases,  they  will  be 
the  only  such  data  available. 
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finally,  a  word  should  bo  said  about  the  applicability  of  the  present  results  to 
real  life.  The  results  obtained  and  the  conclusions  drawn,  are  valid  only  in  so  far 
as  the  Bent  ionospheric  model  is  valid.  While  this  model  seems  to  be  one  of  the 
best  available  for  the  present  purposes,  it  predicts  only  a  monthly  average  varia¬ 
tion  of  the  ionosphere  for  magnetically  quiet  conditions.  Consequently  it  is 
recommended  that  the  conversion  errors  discussed  in  this  report  be  taken  as  lower 
limits  when  being  applied  to  realtime  observations  of  TEC. 
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